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In this communication, we examine a nonlinear model with an impurity emulating a bend. We 
justify the geometric interpretation of the model and connect it with earlier work on models including 
geometric effects. We focus on both the bifurcation and stability analysis of the modes that emerge 
as a function of the strength of the bend angle, but we also examine dynamical effects including 
the scattering of mobile localized modes (discrete breathers) off of such a geometric structure. The 
potential outcomes of such numerical experiments (including transmission, trapping within the bend 
as well as reflection) are highlighted and qualitatively explained. Such models are of interest both 
theoretically in understanding the interplay of breathers with curvature, but also practically in 
simple models of photonic crystals or of bent chains of DNA. 
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I. INTRODUCTION 



In the last decade, intrinsic localized modes (ILMs), or discrete breathers (DBs) as they are also termed, have 
become a topic of intense theoretical and experimental investigation; see e.g., |l| for a number of recent reviews on 
the topic. Per their inherent ability to bottleneck and potentially transport the energy in a coherent fashion, such 
exponentially localized in space and periodic in time entities have have come to be of interest in a variety of contexts. 
These range from nonlinear optics and arrays of waveguides 2] to Bose-Einstein condensates inside optical lattice 
potentials [3| and from prototypical models of nonlinear springs [J| to Josephson junctions |5j and dynamical models 
of the DNA double strand [6|. 

One of the playgrounds that have most recently been added to this long list of ILM applications consists of nonlinear 
photonic crystal waveguides and circuits 0. In connection to this context, an issue that becomes very relevant (see e.g. 
the models developed in 0) is the interplay of nonlinearity and geometryas, typically in photonic crystal waveguide 
arrays, two-dimensional or quasi-one dimensional settings with bends 00 become relevant. 

We should note here that the interplay of nonlinearity and geometry has been increasingly appreciated in the ILM 
literature. From the long range interactions on a fixed curved substrate 0, to lattice-substrate feedback models 
flo| and from lattice junctions with different masses to semi-circular, polymer- like chains an d geometrically 
motivated, bent models of DNA [iH. FLU ITEl IT^| . the geometry can significantly affect the static (inducing e.g., multi- 
stability) and dynamic (causing e.g., a variety of outcomes in the ILM interaction with curvature) properties of the 
relevant lattice model. 

In the present work, motivated by these studies we examine the nonlinear Klein-Gordon variant of a model intro- 
duced recently in the context of the discrete nonlinear Schrodinger equation (DNLS) in [8j|. This is a prototypical 
model emulating the geometry of a lattice bend, by the inclusion in the vicinity of the bend of next-nearest neighbor 
interaction due to the proximity of these neighbors in this context (see e.g., Fig. Q). 

We will examine this bend in the framework of a soft and a hard Klein-Gordon interaction potential, namely the 
Morse and the hard ip 4 potentials respectively. Apart from the inherent interest of the interaction, at the static as 
well as dynamical level, of discrete breathers with the "geometric impurity" induced by the bend, the model may 
be relevant to a number of applications. In particular, the DNLS serves as the envelope wave equation (at the 
discrete level) for Maxwell's equation, hence, for the photonic applications, it may be of more interest to identify 
the properties of the corresponding Klein-Gordon model. Furthermore, this may also serve as a simple dynamical 
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FIG. 1: Schematic presentation of a bend in connection with the discrete equation ©• The parameter 7 stands for the bend- 
induced interaction between the next-nearest neighbors (in the vicinity of the bend), whereas C represents the nearest neighbor 
interaction in the rectilinear chain. 



model for understanding the interaction of a denaturation bubble in DNA [y, [lTj with the local helical geometry of 
the double strand. This is similar in spirit to the earlier study of |l9j| . where the helicoidal geometry was argued to 
induce a non-nearest neighbor interaction across the double strand. Notice, also, that we follow a slightly different 
path than by examining not only static but also dynamic properties of the model and, in particular, the potential 
outcomes of the breather-bend interaction. 

In the next section, we will present the model equation of interest, while in section 3, we will study its linear modes. 
In section 4, we will examine the static properties of the two models, by examining the bifurcations of breathing modes 
as a function of the "bend parameter" 7. In section 5, the corresponding dynamics properties of the breather-bend 
interaction will be presented. Finally, in section 6, we summarize our findings and present our conclusions. 



II. THE MODEL 



In accordance with the above description, the Hamiltonian of the Klein-Gordon chains of interest will be given by: 



H = ^ ( + V ( u n) + 2 C ~ u «) 2 + 7( u «+i - m«-i) 2 <5„,o j , (1) 

thus the bending point is located at the particle n = and implies, as is graphically indicated in Fig. ^ a geometrically 
induced coupling of next-nearest neighbors (nnn) adjacent to the bend site. 

We should clarify here the nature of the configuration and of the dynamics. The displacements u n represent an 
"external" field such as e.g., the electric field in a array of optical waveguides or the stretching of a base pair in a 
DNA chain |g. The chain of Fig. ^should be considered as being geometrically fixed on the plane of the paper. 
The dynamics of the displacements is transverse, i.e., perpendicular to the plane of the bent chain. In such a setting 
the coefficient of the interaction between adjacent sites is determined by the (fixed) geometric proximity of the nodes 
constituting the chain. For a derivation of the discrete model in the first setting and an explanation of the dependence 
of the interaction coefficient on the geometry of the configuration and the inter-nodal distance, see e.g., In this 
context, it is natural to assume that the effect of the bend (due to the geometric proximity of sites n = 1 and n = — 1) 
will be to induce the nnn coupling between them (cf. [lsf). 

V(u) is the on-site potential which is chosen to be of the (soft) Morse type, V(u) = (cxp(— u) — l) 2 /2 or of the 
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(hard) <fi 4 type, V(u) — x 2 /2 + x 4 /A; C is a coupling constant between nearest neighbors (nn) whose interaction will 
be harmonic (i.e., connected by linear springs) for the purposes of the present work. 
From the Hamiltonian (JJJ, the following dynamical equations can be obtained: 

U n + V'(u n ) + C {u n+1 + U n -i - 2u„) + j[(Un - U n - 2 )5n,l + K - U n +2)^rj,-l] = 0. (2) 

A natural bifurcation/continuation parameter that we will use for the purposes of our study is a = j/C. The 
relevant ratio is a natural measure of the relative strength of the different neighbor interactions (the bend-based nnn 
one and the uniformly distributed nn one). This ratio can also be interpreted geometrically (in the sprit of Fig. 
nj. In particular, it can be related to the wedge angle (f> through the relation: a = 1/(2(1 — cos(0))). For the nnn 
approximation to be realistic, ip must be larger than 60°, hence, equivalently, a < 1. Notice that while the geometric 
interpretation of the moment prompts us to typically examine the regime of a £ (0, 1], it is of inherent mathematical 
interest to examine broader parameter ranges, and therefore some of our results below will be presented for values of 
a outside this range. 



III. LINEAR MODES 



Some of the properties of discrete breathers are related to the existence or non-existence of linear localized modes. 
These modes appear as the bend is introduced by modifying the coupling, which is effectively equivalent to the 
introduction of an inhomogeneity in the curvature at the bottom of the inter-site potential pp| |. 

The linear modes can be obtained from the linearized form (around the uniform state u n — 0) of the dynamical 
equations (J5J: 



+ L0 o U n + C{2u n - U n -i - Un+i) + 7[(ltn - U n _ 2 )£„,l + (U n - U n+2 )£ n ,-l] = 0, (3) 



with uj 2 — 1. The linear modes can be calculated using the lattice Green's function (2^, ■ The frequency of the 
localized linear modes (also referred to as impurity modes) is thus given by: 



<„ P - < + Ya = U ° + 2 7 ' (4) 

and the inhomogeneity parameter is related to the frequency of the impurity modes through the relation: 



" = W , (5) 

where the plus sign corresponds to u? mp > w„ + 4C and the minus sign to u>i mp < ui . 

Figure [3 shows the frequencies of the linear modes as a function of a = 7/C and the profile of the impurity modes. 
It can be observed that, for 7 > the impurity mode is above the phonon band and, in consequence, has a zi gzag 
vibrational pattern (i.e., a staggered mode). If 7 < 0, the mode is below the band and its sites oscillate in phase [2(1 . 
It is worth remarking that the impurity modes are anti-symmetric modes whose central lattice site is at rest. 



IV. STATIC RESULTS: STATIONARY BREATHERS AND THEIR BIFURCATIONS 



Discrete breathers can be calculated using well-known techniques based on the concept of continuation from the 
anti-continuous limit |23| . Furthermore, since the bending acts as an inhomogeneity, the system has lost its integer 
shift invariance and the properties of the breathers critically depend on the site where their center is located (at 7 = 0). 
This fact leads to the existence of saddle-node bifurcations through which some of the solutions can disappear. Similar 
bifurcations have been previously observed in bent chains of oscillators [l3l Il4j. This has also been the main focus, 
for DNLS type settings, of 

These bifurcations need a relatively high value of the coupling (sufficient cross-talk between the neighboring sites) 
to be observed. Otherwise, breather solutions will exist for very large intervals of 7 in every site of the lattice. Hence 
we restrict ourselves to C"s of the order of (typically) 0.1 to render these bifurcations tractable. 

We now proceed to examine the bifurcation diagrams (as a function of a) for both soft and hard on-site potentials 
and a harmonic inter-site potential. 



4 




FIG. 2: (Top) Linear modes spectrum with respect to the bending parameter a = 
(Bottom) Profile of the impurity mode for a > (left) or a < (right) and |a| = 2.5. 
modes is non-negligible only for high values of a. 



- y/C for a coupling constant C — 0.20. 
Note that the localization of the impurity 



A. Morse potential 

Figure |3| shows the bifurcation diagrams for the case of a Morse on-site potential and couplings C = 0.13 and 
C = 0.26 and frequency lj-^ = 0.8, and Figure 0] shows the spatial profile of the solutions corresponding to a number 
of branches for a = 0.0005. 

For C = 0.13, the branch corresponding to the breathers centered at n — merges with the breathers originally 
centered at n = 0.5 at a ps 0.00115. The branches centered at n = 1.5 and n = 2 disappear through a saddle-node 
bifurcation at a ~ 0.00057. This phenomenon is also observed for the branches centered at n — 2.5 and n — 3 at 
q ps 0.00674. For a < the annihilations are observed between different branches of solutions. In particular, the 
branches centered at n = 0.5 and n — 1 annihilate at a w —0.00049 and the ones centered at n = 2 and n = 2.5 
cease to exist through the saddle-node bifurcation occurring ata« —0.00120. It should, however, be noted that some 
branches (see e.g., the branch centered at n — 1 for a > or the one centered n = — 1.5 for a < 0) never annihilate. It 
is also interesting to note that all of these branches of solutions are unstable, as both site-centered and bond-centered 
breathers are unstable in the straight chain (for this value of the nn coupling). However, the branch resulting when 
the centered at n = and n — 0.5 merge is stable. 

The branch centered at n = 1.5 becomes stable for a < 0.002. This change of stability has its origin in a pair of 
Floquet multipliers that collide in the unit circle at A = 1. Notice that this behavior, unlike the typical branches 
discussed above does not involve a saddle-node bifurcation (or for that matter an exchange of stability with another 
static breather branch). 

For C = 0.26 the saddle-node bifurcations occur between different solutions from the case considered above. For 
a > 0, the branches centered at n = 0.5 and n — 1 collide and disappear at a w 0.00072; the branches centered 
at n — 2 and n = 2.5 cease to exist at a w 0.00106 while the ones centered at n — 3 and n = 3.5 terminate at 
a w 0.00342. The branch centered at n = merges with the branch centered at n — 0.5 for a < and in particular 
for a « -0.00187. 

In addition, the branches centered at n = 1.5 and n = 2 annihilate at a ps —0.00204 and the ones centered at 
n — 2.5 and n = 3 collide at a ps —0.00157. It is interesting to observe that similarly to what was found in [8|, in 
a different setting, there is an asymmetric mode that persists for a < 0, whose energy is lower than the symmetric 
bend mode. Hence, for a < in this case, we observe a symmetry breaking effect that leads to an asymmetric ground 
state of the system. Notice that similar asymmetric modes can be found (and identified to be potentially stable) in 
continuum models with localized impurities; see e.g., 0. 

It should be noted here that for this larger value of the coupling, the branches corresponding to site centered 
solutions are stable for an infinite lattice, whereas the bond-centered solutions are unstable for the rectilinear chain. 
This stability is inherited by the branches of the bent chain. Motivated by the change of stability at the "un-bent" 
limit (7 = 0), occurring as a function of C (e.g., notice above the different stability of this limit for C = 0.13 and 
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FIG. 3: Bifurcation diagram for the solutions with a Morse on-site potential. The energy of the breathers is plotted as a 
function of the bending parameter a = j/C. The coupling constant is C = 0.13 (top) or C = 0.26 (bottom). The numbers 
indicate the site where the solutions are centered; an integer number corresponds to a site-centered solution, whereas a half- 
integer number corresponds to a bond-centered solution. Stable solutions (for an infinite system) are represented by full lines 
whereas unstable solutions are represented by dashed lines. 



C = 0.26), we briefly study the stability of the rectilinear chain limit, as a function of C. To illustrate the stability of 
the 1-site and 2-site modes in the rectilinear chain (discussed in part in |l6j]). we examine their Floquet multipliers, 
both by means of the argument angles as well as giving their absolute values in Fig. [5] Notice, in particular, for 
the 1-site modes that there are two inversions of stability, namely one at C w 0.1297 and another for C « 0.2382. 
For C > 0.2382, the 1-site modes are stable but with interspersed size-dependent windows of instability (see e.g., 

miMm). 

Finally, let us mention that the branch that originates by the merging of the n = and n = 0.5 branches is unstable. 
Figure |H1 shows an example of the spatial profile and the Floquet eigenvalues of an unstable state for a < 0, namely 
a breather (originally) centered around n = 1. Notice however that the instability in this case is a result of the finite 
size of the lattice. 
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FIG. 4: Spatial profiles of the breathers centered at n = (a), n = 0.5 (b), n = 1 (c) and n = 1.5 (d) for a — 0.0005. The 
results correspond to a chain of oscillators with an on-site Morse potential with parameters Wb = 0.8 and C = 0.26. 




FIG. 5: Evolution with respect to the coupling constant of the Floquet multiplier arguments 9 (left) and the modulus of the 
Floquet multipliers (right) for a 1-site (top) and a 2-site (bottom) breather with Morse on-site potential with frequency = 0.8. 
The 1-site breather is stable for C € (0,0.1297) and for C S (0.2382,0.4) it recovers its stability except for size-dependent 
instability bubbles via oscillatory and subharmonic bifurcations. The 2-site breather is stable for C £ (0.1300, 0.2381). 



B. cj> 4 potential 

The most significant difference of the cf> 4 case with respect to the Morse one is that in the former (i.e., for a hard 
potential), the breathers bifurcate from the top of the continuous spectrum and as a result the corresponding spatial 
profiles are staggered 0. Hence, in this case the main solutions of interest are 1-site or 2-site breathers with the 
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FIG. 6: (a) Spatial profiles of the breathers centered at n — 1 for a = —0.003 and the same characteristics to those of Figure 
^] (a) shows the spatial profile of unstable solution; (b) represents the Floquet multipliers corresponding to this state. The 
instabilities that can be appreciated near 8 = n are due to the collision of extended eigenmodes and disappear in the case of 
an infinite lattice. 
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FIG. 7: Bifurcation diagram for the solutions with a (f> 4 on-site potential. The energy of the breathers is plotted as a function 
the bending parameter a — 7/C. The numbers indicate the site where the solutions are centered; the half-integer numbers 
actually correspond to 2-site breathers, but the notation has been kept in consonance with Figure Stable solutions are 
represented by full lines whereas unstable solutions are represented by dashed lines. 



adjacent sites oscillating in anti-phase. 

The hard c/> 4 potential has been considered with a coupling constant C = 0.09 and solutions with frequency oj^ = 1.2. 
For these parameter values (notice that C is comparable to the values used in Fig. |3J), the saddle- node bifurcations 
of interest occur at considerably larger values of a for the present model in comparison to the Morse case; see e.g., 
Fig. In particular, the branch of the on-site breathers centered at n = merges with the inter-site branch centered 
at n = 0.5 at a w 0.516. Similarly, the on-site solutions centered at n = 2 and n = 3 collide (and disappear) with 
the corresponding inter-site ones centered at n = 1.5 and n — 2.5 at a w 0.258 and a m 0.590, respectively. However, 
in this case, the branch centered at n = 1, becomes the asymmetric, lowest energy, stable ground state that persists 
for any value of a. For a < 0, the on-site breathers centered at n = 1 annihilate with the inter-site ones centered at 
n = 0.5 at a w —0.255. Notice that here also on-site solutions are stable, while inter-site ones are unstable for 7 = 0, 
but the branch emerging from the n = and n = 0.5 is unstable. 

Some of the breather solutions for a = 0.1 are shown in Figure [S] and the profile of the ground state for a > 
together with its Floquet eigenvalues in Figure 
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FIG. 8: Spatial profiles of the breathers centered at n — (a), n = 0.5 (b), n = 1 (c) and n = 1.5 (d) for a — 0.1. The 
results correspond to a chain of oscillators with an on-site (f) 4 potential with breather frequency = 1-2 and coupling constant 
C = 0.09. 
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FIG. 9: (a) Spatial profiles of the breathers centered at n = 1 for a = 0.5 and the same characteristics to those of Figure [S] 
This solution corresponds to the stable ground state, (b) represents the Floquet eigenvalues of this solution. 

C. Instabilities: Switching and Breather Mobility 

As explained above, the ground state is not always the breather centered at n = 0. In particular, for a Morse 
potential, C = 0.26, Wb = 0.8 and \a\ large enough, the ground state for a < is the breather centered at n = 1 (at 
least in the infinite domain limit, since in the finite domain case, it may be unstable as shown in Fig. On the other 
hand, for a > the ground state is still the breather at n = 0. Thus, if an unstable solution (say, the n = 0-centered 
for a < or the n = 1.5 for a > 0) are perturbed, they are likely to switch to the ground state, in a similar fashion to 
the phenomenon observed in [l^. However, as Figure ITU1 shows, this may not always be the case. In particular, the 
breather oscillates between the sites n = and n = 3.5 (i.e., in the vicinity of the stable ground state) for a < 0, while 
for a > 0, it is set in motion. A possible explanation relies on the fact that the eigenmode leading to the instability 
is actually a pinning mode, that is responsible for breather mobility (see also Section [Vj. Hence a perturbation with 
the appropriate sign (i.e., causing the breather to move away from the ground state of n = rather than toward it) 
along this eigendirection can lead to breather motion. 

For a <fi 4 potential, C — 0.09, cub = 1.2 and |a| large enough, the ground state for a > is the breather centered at 
n=l, while for a < 0, it is not well-defined. 
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FIG. 10: Time evolution of an unstable breather slightly perturbed for a > (a) and a < (b). 
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FIG. 11: Time evolution of an unstable breather slightly perturbed for a > 0. 



Figure ITT1 shows the evolution of a 2-site breather centered at n = 2 with a = 0.5. It is observed that the breather 
can switch to the 1-site breather centered at n — 1, i.e., the ground state. In this case, no moving localized excitations 
arise, but rather the previously suggested switch occurs. This may be attributed to the fact that for the hard (f> 4 
potential, there do not exist pinning modes that could potentially lead to breather mobility (see also [28j). 



V. DYNAMIC RESULTS: MOVING BREATHERS AND BREATHER-BEND INTERACTIONS 

In this section we examine how the local geometry of the bend can influence the motion of a breather in its vicinity. 
In particular, the moving breathers are launched towards the bend vertex following the marginal mode method 28] . 
The latter consists of adding to the velocities of the static breather a perturbation of magnitude A, collinear to the 
direction of a specific linear localized mode, namely the eigenmode corresponding to translation. Then, the kinetic 
(translational) energy of the breather can be defined as K = A 2 /2. As Ref. [l5| shows, the moving breather can 
be described as a quasi-particle of mass m* defined through the relation K = A 2 /2 = m*v 2 /2, with v being the 
translational velocity of the breather. The mass is a measure of the inertia of the breather to external forces. 

Another useful concept for our study is the breather energy center, which is a measure of the position of the localized 
structure. It is defined through the relation X = (52 n ne n )/E, i.e., as the center of mass with respect to e„, the 
energy density at the n-th site. E = ^ n e n is the total energy of the breather [Tfij . 

The local "coupling inhomogeneity" induced by the bend, should result in features similar to those observed in Ref. 
[2fj|. In particular, for a soft potential, if a > 0, the bending should act, in principle (see also below), as a potential 
barrier, i.e., the breather can cross the vertex as long as its translational energy is above a threshold value. If this 
condition is satisfied, the breather decelerates when crossing the bend, but eventually recovers its initial velocity. If, 
on the other hand, a < 0, then the corresponding behavior is more complex: for small values of \a\, the bend acts as 
a potential well, accelerating a moving breather approaching it. However, for larger values of the parameter, trapping 
and even reflection become possible. For a hard potential, such as 4 , the behaviors for a > and a < are reversed 

We will focus on the a > case and consider, as before, a chain of oscillators with 1) the Morse on-site potential, 
as well as one with 2) the </> 4 on-site and inter-site potential. Thus, for the case of the Morse (soft) potential, the 
bending will act as a potential barrier, while for the </> 4 (hard) potential, the well/trapping/reflection regimes will be 
examined. 
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FIG. 12: Evolution of the energy center (left) and the moving breather (right) for a Morse on-site potential with parameters 
uJb = 0.8 and C = 0.13. Figures (a) and (b) correspond to a reflection (a = 0.008 and K — 0.00174), (c) and (d) to a refraction 
(a = 0.008 and K = 0.00180) and (e) and (f) to a refraction through the site n = — 1 and a reflection at n = 1 leading to a 
trapping (a = 0.008 and K = 0.00392). 



A. Morse Chain: Bend-Induced Potential Barrier 

As mentioned earlier, for a soft on-site potential and a > 0, the bending acts as a potential barrier, i.e., if the 
translational energy of the breather K is above a threshold value U c , the breather crosses the bending. Otherwise, 
the breather is reflected. Figure IT21 shows the time evolution of a reflected (top panel), a refracted (middle panel) and 
that of a trapped breather (bottom panel). In the top panel (for K < U c ), the breather is immediately reflected from 
the potential barrier, emerging with (approximately) the same speed, but propagating in the opposite direction. In 
the middle panel of the figure, it can be observed that, for K > U c the breather decelerates when it approaches the 
n = — 1 site, accelerates when it crosses this point until it reaches the bend vertex (n = site); it decelerates again 
when approaching the n = 1 site, accelerating when it crosses this point. 

Figure^]shows the critical energy U c necessary to cross the bend. An important point to indicate is the existence of 
two distinct values of the critical energy U c ± , that correspond respectively to the kinetic energy threshold for crossing 
the site n = ±1. When U c - < U c + holds, for intermediate values of U- < K < U+, the loss of energy due to phonon 
radiation in the breather-bend interaction, disallows the moving breather from crossing the bend. Hence, if the two 
thresholds are different, the breather can be trapped oscillating between the n — —1 and n = 1 sites; an example is 
shown in the bottom panel of Fig. As is naturally expected, the energy loss decreases when a is small (i.e., as we 
approach the rectilinear, integer-shift invariant chain), and U c - = U c + is practically fulfilled in this parameter range, 
e.g., for a < 0.01 in Fig. U3 

The potential barrier U — U(X) can be calculated using the procedure described in ^5|: if K < U c , the value of 
K is fixed and X Q is the point where the breather is reflected (turning point), then K = U(X ). 
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FIG. 13: Dependence of the critical energy U c on a for a Morse on-site potential with parameters uib = 0.8 and C = 0.13. Note 
that for sufficiently large a, the curve separates in two parts corresponding to U c + and U c - (see also the relevant discussion in 
the text). 



If, on the other hand, K > U c , the barrier can be calculated through the formula: 



where v is the initial translational velocity of the breather and v(X) is its (numerically computed) velocity at the 
point X. If the mass of the breather is constant, the barrier calculated for K < U c and K > U c should coincide, 
except for the oscillations that appear in the second case with respect to the first one. The origin of these oscillations 
relies on the non-uniformity of the instantaneous translational velocity of the breather due to the discreteness of the 
system. Thus, to obtain the actual shape of the barrier, the turning point approach should be used. 

However, this method can only be applied in our case for X < — 1 (and hence by symmetry for X > 1) as for K < U c , 
the breather cannot cross over to the region — 1 < X < 1. In order to overcome this drawback, we have approximated 
the potential barrier created by a single inhomogeneity with a Gaussian shape: U(X) = aexp(— b(X — X n ) 2 ), where 
X Q is the point where the maximum of the barrier is located (which need not be the inhomogeneity point) |8L I26L l30| . 
In our bend setting, there are two points where the equation is inhomogeneous with respect to the rest of the chain. 
Hence, it is natural to approximate the actual barrier with the linear superposition of two Gaussian barriers, i.e., 



U(X) = aexp{-b(X - X^) 2 ) +exp(-b(X - Xrf), x ±1 =±{l + S), (7) 

with 6 being the difference between the turning point at K ss U c and the "impurity site"; e.g. for a breather launched 
from the left of the bending, 6 = X(U C ) + 1. Thus, equation Q becomes: 



U(X) = 2aexp{-bX 2 ) exp(-6(l + S) 2 ) cosh(2&X(l + 5)). (8) 

Figure ITU shows the barrier calculated using equations © and JSJ. Parameters a, b and 6 are chosen through the 
Gaussian fit of the barrier points for \X\ > 1. Comparing the fitted profile of Eq. (JSJ with the one obtained for a 
supercritical case of K > U c indicates very good agreement, which in turn supports the assumption of Eq. Q. 

Notice, however, that the such a restricted one-degree-of-frcedom picture cannot capture phenomena such as emis- 
sion of radiative small amplitude waves during the crossing of the bend, and consequently cannot display a phe- 
nomenology as rich as the original problem (e.g., it cannot exhibit the trapping type phenomena discussed above). 
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FIG. 14: Potential barrier experienced by the moving breather in the chain with a Morse on-site potential, when it reaches 
the bending point. Parameters of the breather are uib = 0.8, C = 0.13 and 7 = 0.08. The thick line corresponds to the fit to 
equation JHJ and the slim line to equation © for K — 0.0018. The oscillations of the last curve are due to the non-uniformity 
of the translational velocity of the breather. Parameters of the fitting to equation JSJ are: a — 0.001768, b — 2.7068 and 
8 = 0.08754. 



B. 4> 4 Chain: Well/Trapping/Reflection Regimes 

For the case of a hard potential with a > three different behaviors can be observed when a increases (see Figure 

1. For small values of a, the breather accelerates when reaching the bend points, and eventually emerges from the 
bend recovering a constant velocity (but not the same as the initial one due to radiative emissions). That is, 
the bending acts as a potential well. 

2. For larger values of a, the breather is trapped at the n = — 1 site. 

3. For even larger values of a, the breather is reflected and the particle at the vertex remains excited. This 
phenomenon can be also called partial trapping, since a fraction of the original breather energy remains localized 
at the bend. 

4. Finally, beyond a third critical point for a, the breather is reflected from the bend, and, in a different fashion 
to the barrier case, the breather is not refracted for any velocity. 

It should be noted that the boundaries between partial trapping to reflection and well to trapping regimes are not 
well-defined. In the first case, the amplitude of the vertex excitation decreases when \a\ increases. Due to this fact, 
it is difficult to establish the borderline between those regimes. In the second case, the breather is trapped if its 
translational energy is below a threshold value. However, for energies above this value, the bend acts as a potential 
well. Thus, the trapping takes place without need of energy loss when the breather crosses the n = — 1 site, that is, 
it occurs as long as the translational energy is below the threshold commented above. These facts are illustrated in 
Figure ED 

VI. CONCLUSIONS 

In this communication, we have introduced a discrete Klein-Gordon chain model variant, emulating the existence 
of a geometric bend at a lattice site. 

We have examined the static properties of the localized in space, time periodic solutions that exist in this setting, 
starting from the limit where the chain is rectilinear. In the latter case, from the integer shift invariance and 
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FIG. 15: Evolution of the energy center (left) and the moving breather (right) for 4 on-site and interaction potentials with 
parameters uj^ = 3 and C = 0.5352. Figures (a) and (b) correspond to a refraction (a — 0.2), (c) and (d) to a trapping 
(a — 0.5) and (e) and (f) to a reflection (a — 4). 




numerous earlier studies [l|, we know that site-centered and bond-centered solutions will exist and their stability 
can be established. Hence, using continuation from the "straight" limit into the "bent" regime, we are able to trace 
the branches of corresponding solutions and identify their saddle-node bifurcations leading to the termination of the 
corresponding branches. However, we have also identified symmetry breaking effects in which one of the asymmetric 
branches does not collide with one of its neighboring counterparts to terminate its existence, but rather survives 
throughout the bend-parameter strength continuation, and becomes the lowest energy state of the lattice. Excitation 
of unstable bend modes in such cases will lead to trapping around (i.e., switching to) such an asymmetric state, but 
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in the presence of pinning modes, moving breathers may result. Hence, we conclude that sufficiently strong bends 
may favor an asymmetric localization of energy, with respect to the bend center. 

We have also studied the dynamic properties of the interaction of such bends with moving localized modes, when 
the latter are scattered off of a bend. We have found that for soft potentials, the bend operates as a potential barrier 
allowing transmission for supercritical and reflection for subcritical values of the localized excitation's initial speed. 
However, there can also be narrow intervals of trapping (not captured by the potential barrier picture). On the other 
hand, for hard potentials, a more complex and partly initial condition (i.e., kinetic energy) dependent picture emerges. 
For small bend strengths, the inhomogeneity acts as a potential well, but for higher ones it can lead to trapping and 
even to reflection. 

In this numerical study, we have mainly aimed at presenting the relevant phenomenology, classifying it for Klein- 
Gordon chains with different potentials, and providing a qualitative explanation of the relevant findings. It would 
naturally be of interest to attempt to obtain more quantitative, theoretically predicted, estimates for the transitions 
observed herein. A variational viewpoint may prove to be useful in this context. 
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